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Abstract 

We give upper bounds for the number of spin ^ particles that can be 
bound to a nucleus of charge Z in the presence of a magnetic field B, in- 
cluding the spin-field coupling. We use Lieb's strategy, which is known to 
yield N c < 2Z + 1 for magnetic fields that go to zero at infinity, ignoring 
the spin-field interaction. For particles with fermionic statistics in a homo- 
geneous magnetic field our upper bound has an additional term of order 
Z x min{(5/Z 3 ) 2 / 5 ,l+ [lnCB/Z 3 )! 2 }. 



1 Introduction and main result 

Let Hn,z,a be the Hamiltonian for N identical particles with spin 5 in the 
Coulomb field of a nucleus of charge Z and in a magnetic field B = curl A, 

i=l v 1 %w i<j 1 Jl 

with 

h a = [<*j ■ HVj + A( Xj ))] 2 = (-iVj + A( Xj )) 2 + B • o> (2) 

Here A 6 /^(IR 3 ; IR 3 ) is the magnetic potential and cr = (a x , a y , a z ) are the 
usual Pauli matrices. The Hamiltonian (0) acts on the fermionic (respectively 
bosonic) subspace of (gT £ 2 (R 3 , rfx; C 2 ). We will assume that the ground state 
energy 

E(N, Z, B) = inf spec H N ^ A (3) 

is finite. (Note that the energy depends only on B because of gauge invari- 
ance.) A sufficient condition for this is B G £ 3 ^ 2 + C°°, because is this case 
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|B| + Z/\x.\ is relatively bounded with respect to —A, and inf spec Hn,z,a > 
inf spec J2i=i ~ |B(xj)| — Z/|xj|) by the diamagnetic inequality More- 
over, we will only consider magnetic potentials A such that the energy E(N, Z, B) 
is monotonously decreasing in N for fixed Z. This is in particular the case for a 
homogeneous magnetic field. We are interested in the maximal number of parti- 
cles that can be bound, i.e. the largest N such that E(N, Z, B) is an eigenvalue. 
We will denote this "critical" N by N c , suppressing the dependence on Z and B. 
For simplicity we will restrict ourselves to considering identical particles. 
Alternatively, one could define the critical particle number by 

N c = max {N\E(N, Z, B) < E(N — 1, Z, B)} . (4) 

With this definition E(N, Z, B) is certainly an eigenvalue if N < N c , so N c < N c . 
Hence any upper bound to N c is also an upper bound to N c . 

It is well known that magnetic fields, at least homogeneous ones, enhance 
binding. In H it is shown that every once negatively charged ion (i.e. N = 
Z + 1) has an infinite number of bound states in the presence of a homoge- 
neous magnetic field of arbitrary field strength B. And in |3| the lower bound 
liminf^_ >00 (A c /Z) > 2 as long as B/Z 3 — > oo is given, which is in contrast to 
asymptotic neutrality in the absence of magnetic fields |5|]. 

We will use Lieb's strategy to derive an upper bound on N c . The difference 
between our considerations and |T[ is the coupling of the spin to the magnetic 
field, i.e. the term B • a in the Hamiltonian. Without this term, Lieb derived the 
bound N c < 2Z + 1 for any bounded A that goes to zero at infinity. 

Our result is as follows: 

THEOREM 1 (Upper bound on N c ). Under the conditions stated above, 

1 E(N C , Z, B) - E(N C , kZ, B) 

for all values of k > 1 . 

Note that since the ground state energy is superadditive in N, E(N, Z, B)/A 
is bounded by some some function independent of N. Moreover, the best bound 
in (0) is achieved in the limit k \ 1, which exists by concavity of E(N, Z, B) in 
Z. 

To apply Theorem [L] to the case of fermionic electrons in a homogeneous 
magnetic field B = (0,0,5), one needs upper and lower bounds to the ground 
state energy. These were derived in and and are given in section ^4.1| . The 
result is the following: 

THEOREM 2 (Maximal ionization for fermions). Let H N z ,a be the re- 
striction of (ji[) to the fermionic subspace, and let B = (0, 0, B). Then, for some 
constants C\ and Ci, and for all values of B > and Z > 0, 

N C <2Z + 1 + dZ 1 ^ + C 2 Z min {(5/Z 3 ) 2/5 , 1 + I ln{B/Z 3 )\ 2 } . (6) 
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Of course we do not believe that these bounds are optimal. One might assume 
that Lieb's bound N c < 2Z+1 holds also in this case, at least for large Z (compare 
with the lower bound in || stated above), but it remains an open problem to 
show this. However, Theorem |2| improves a result obtained in ||, which states 
that N c < 2Z + 1 + cB x l 2 for the Hamiltonian (|1|) restricted to some special wave 
functions in the lowest Landau band, which reduces the problem to an essentially 
one-dimensional one. 

One might ask how the Pauli principle affects the result in Theorem 0. It 
turns out that the analogue for bosonic particles is the following: 

THEOREM 3 (Maximal ionization for bosons). Let H nz ,a be the restric- 
tion of to the bosonic subspace, and let B = (0, 0, B). Then for some constant 
C3 and for all B > and Z > 



Z 

N c < 2Z + 1 + — min 




in, | 



(7) 



In the next section we will give the proof of Theorem [I]. In section |3| several 
possible generalizations are stated, and in section |] the necessary energy bounds 
for the case of a homogeneous magnetic field are given, which will prove Theorems 
| and |. 

2 Proof of Theorem J] 

Let $ be a normalized ground state for Hn,z,a- Assume, for the moment, that 
(*||xjv|^) is finite. Then 



E(N,Z,B)(\x N \V\V) = (\x N \*\H NtZ>A V) 

4 I y , ■ I - — ' I v ■ — v > ■ I ; 

(8) 



Xiv Xj - Xiv 



Because r = J \I/*\I/|xjv|c/xjv is an acceptable trial density matrix, we can use the 
variational principle to conclude that 

(\x N \*\H N - ltZtA \V) > E(N - 1, Z, B)(\x N \*\V). (9) 

By assumption the energy is monotonously decreasing in N, so 

< (E(N, Z, B) - E(N - 1, Z, B)) (|xjv|tf |tf ) < 0. (10) 
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Now using the demanded symmetry of we get 



X/v| 



|Xjy| + |Xj 
|Xj - Xjv| 



(the strict inequality follows from the fact that {(x, y), |x — y| = |x| + |y|} has 
measure zero), so (TTOh gives 



z> 1 -(N-i) + (\x N \m i A¥)- 



(12) 



As in [0] we have for any positive function y?(xjv) 



(<p^f\H^\^) = & 1/2 M ( H 



(TV) 



2<p 



-zRe (^/ 2 ^|^ ■ i-iV + A)^ 2 V). 

■ <p 



(13) 



Now (\x.n\^?\H^^>) is certainly real, because all the other quantities in equation 
(PP are real. So choosing y(xjv) = |x/v| in equation (TT3]) we get 



(\* N \m!P*) = (\xN\ 1/a *\(HP-^ 

V ^|Xy 

Using that > equation (|12"D reads 



1/2 *). 



(14) 



/V< 2Z+l + -(^||x Jv r 1 *}. 



(15) 



Moreover, 



(^Hxjvl-^Xfe-l) 



< 



NZ 
1 



iV,Z,A _ Hn kZ,A 



(E(N, Z, B) — E(N, kZ,B)) 



(16) 



so we arrive at the desired bound for N c . 

Throughout, we have assumed that (|xjv|\l/|\l/) is finite. A priori, this need 
not be the case. However, one could arrive at the same conclusions using the 
bounded function <£> £ (xat) = |xjy|(l + £|xat|)~ 1 instead of |xjv| in @, and letting 



e — >• at the end (see [0]). Note that 

2 



Vy? e 



< 



|x| 2 (l + £:|x|) 2 |x|yj e (x) 



;i7) 



so our conclusions remain valid. 
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Remark 1. Instead of ignoring the kinetic energy in (|14| ) one could use the 
operator inequality 

HV + Af-^>0 (18) 

to conclude that 

N c <2Z + l-2(*\\x N \B(x N )-(r N *). (19) 

This may especially be of interest if |B(x)| < &|x| _1 for some constant b. And 
for B = Lieb's bound N c < 2Z + 1 is reproduced. 

3 Generalizations of Theorem [I] 

As in [0] several generalizations of Theorem [I] are possible: 

• One can allow different statistics than the bosonic or fermionic one, or even 
consider independent particles. Moreover, the particles could have different 
masses and charges. 

• Hartree- and Hartree-Fock theories can be treated in the same manner. 

• One can replace the Coulomb interaction (everywhere) by some positive 
v(x) = l/w(x), with w satisfying 

w{x-y)<w{x)+w{y), (20) 

and for some constant C 

|Vuf < C. (21) 

Looking at the proof of Theorem |l| we see that these two properties are 
really what we needed. 



4 Application to homogeneous fields 

We will now apply Theorem [l] to the case of a homogeneous magnetic field B = 
(0, 0, B), and prove Theorems and [5]. The magnetic potential, in the symmetric 
gauge, is given by A = ~B x x. The energy in this case will be denoted by 
E(N, Z, B) = E(N, Z, B). To derive explicit bounds on N c , we need upper and 
lower bounds to the ground state energy of ([!]). However, since we are not trying 
to give the optimal constants, the upper bound E(N, Z, B) < will suffice for 
our purposes. So we will concentrate on the lower bounds. We will distinguish 
between the fermionic and the bosonic case. Throughout, every fixed constant 
will be denoted by C, although the various constants may be different. 
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4.1 The fermionic case 



In || and || the following lower bounds on the ground state energy of (]!]) were 
derived: 

LEMMA 1 (Lower bounds on the fermionic energy). Let A = N/Z. The 
ground state energy of ([Jj restricted to the fermionic subspace satisfies: 



(a) IfB< CZ 4 ' 3 then 

E(N, Z, B) > -CZ 7/3 X 1/3 (1 + CX 2/3 ) . 

(b) IfB> CZ 4 ' 3 then 

E(N, Z, B) > -CZ 9/5 X 3/5 B 2/5 (1 + C\- 2/5 ) . 

(c) IfB> CZ 2 then 



E(N, Z, B) > -CNZ 2 1 



In 



A 1 / 2 V Z 3 



(22) 



(23) 



(24) 



Remark 2. Part (c) of Theorem || follows from omitting the repulsion terms in 
Theorem 1.2 ("Confinement to the lowest Landau band") in and then using 
the bound (4.11) there. Although this theorem is applicable for B > CZ 4 / 3 , with 
an additional error term, the result is simpler for B > CZ 2 . However, the bound 
( |24| ) is only of interest for B > CZ 3 , because for smaller B (^3) is more useful. 

Using (HD and the bounds in the preceding Lemma we find that A c = N c /Z 
satisfies 



A c 
A, 



< 2 

< 2 



z- 1 
z- 1 



C\fl 3 Z- 2 ' 3 (l + C\ 2 J 3 ) if B<CZ 4 ' 3 , 
C \-^ (^-\ (l + CA c - 2/5 ) if B> CZ 4 ' 3 , 



A c < 2 + Z- 1 + C 1 



In 



C (B 



- 



\l' 2 \z 3 



1/2 



+ 1 



if B > CZ' 



(25) 



Putting together these three bounds, we obtain the result stated in Theorem ^. 
Note that these bounds imply in particular 



limsup— <2 if B/Z 3 

Z^oo Z 



0. 



(26) 
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4.2 The bosonic case 

To get a lower bound the bosonic energy we will first omit the positive repulsion 
terms in ([!]). By scaling the variables Xj — > Z~ x Xi we see that 

E(N, Z, B) > NZ 2 e(B/Z 2 ), (27) 

where e(b) is the ground state energy of hydrogen in a homogeneous magnetic 
field of strength b. For small b, we will use the diamagnetic inequality, which 
implies 

e(&) >_!-&. (28) 

A large b expansion of e(b) is given in ||. From there we get the following lower 
bound: 

LEMMA 2 (Lower bound for the hydrogen energy). For large enough b 
the ground state energy of hydrogen, e(b), satisfies 

e(6)>-i(ln6) 2 (l + ^. (29) 

Note that the bosonic energy is at least of order NZ 2 , even for small B. 
Therefore the contribution to (|5|) is always at least O(Z), in contrast to fermions, 
where the energy is of order N l ^Z 2 for small B (this is the reason for the addi- 
tional factor Z 1 / 3 in ffl)). Setting k = 2 we arrive at the bound given in Theorem 

i 

We remark that since the bosonic energy is always less than the fermionic 
energy, Theorem |3] holds also for fermions; but the bound stated there is certainly 
worse than the one given in Theorem |2|. 
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